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Abstract 

Simple derivation of the condition for the transition point from absolute insta- 
bility of plane dark solitons to their convective instability is suggested. It is 
shown that unstable wave packet expands with velocity equal to the minimal 
group velocity of the disturbance waves propagating along a dark soliton. The 
growth rate of the length of dark solitons generated by the flow of Bose-Einstein 
condensate past an obstacle is estimated. Analytical theory is confirmed by the 
results of numerical simulations. 
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1. Introduction 

Instability of dark solitons with respect to transverse perturbations is well 
studied both theoretically and experimentally (see, e.g., review articles [T1IHI3] 
and references therein). Qualitatively, it is caused by decrease of the energy 
per unit length of a soliton with increase of the local velocity of the disturbed 
element. As a result, the stretched segments of the soliton acquire increased 
velocity and hence the soliton's local curvature grows with time which leads 
to breaking of the soliton followed by formation of vortex pairs. However, this 
mechanism of the dark soliton instability was studied mainly when the dark 
soliton propagates in a quiescent medium. For example, this situation is realized 
in dynamics of Bose-Einstein condensate (BEC), when dark solitons are formed 
by means of density or phase engineering. However, the situation can change 
drastically, if there is fast enough flow of the condensate along the dark soliton. 
In particular, as was found in [1], just this happens when dark solitons are 
generated by a flow of BEC past an obstacle. In this case dynamics of BEC 
is described (in standard non-dimensional units) by the Gross-Pitaevskii (CP) 
equation 

iV-t + - (|^|2 - 1)^ = C/(r)V', (1) 
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where ip{r,t) is the condensate wave function and U{r) the potential of the ob- 
stacle. This equation was solved numerically in ^ with the boundary condition 



ip\s ~ cxp(iMa;) 



(2) 



at the boundary S far enough from the obstacle, and it was found that there is 
a critical value of the Mach number of the incident flow, 



so that for M < Mcr the vortices are generated downstream the obstacle whereas 
for M > Mcr the dark solitons are stretched from the shadow behind the ob- 
stacle. This difference in the flow behavior was explained in [S] as a transition 
from absolute instability of dark solitons to their convective instability, so that 
for M > Mcr the unstable disturbances of the dark soliton are convected by the 
flow along it with such a velocity that they cannot destroy the soliton at any 
finite distance from the obstacle after long enough time and, hence, the soliton 
increases its length with time. In this letter we shall introduce the notion of 
velocity of propagation of the front of the unstable wave packet into undisturbed 
region of the dark soliton. This will permit us to give simple physical interpreta- 
tion of the condition for the absolute/convective instability transition found in 
[5] and to determine quantitatively the growth rate of oblique solitons length in 
the experiments with flow of condensate past an obstacle what is quite topical 
in connection with recent experimental observation [B] of oblique solitons in the 
flow of the polariton condensate. 

2. Absolute and convective instability of dark solitons 

To simplify the notation, equation ([T]) is written in such units that the back- 
ground uniform density is equal to unity and unit of velocity is equal to the 
sound speed (c = 1). Let the dark soliton propagate along x axis with velocity 
V , then the corresponding dark soliton solution of ([T]) (with U = Q) has the 



ip{x,y,t) = i:,{x-Vt) = Vl - tanh[Vl -V^ix-^ Vt)] + iV. (4) 



Mcr = 1.43, 



(3) 



form 



If we represent the condensate wave function as 



i'{x,y,t) = V "(a^: t) exp{i(t){x, y, t)) 



(5) 



then the condensate density and the flow velocity. 



(6) 



are given for the soliton solution by the formulae 



n = ns(x — Vt) = 1 



cosh^[Vl-V^{x -Vt)] 



(7) 
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Stability of such a solution was studied in [7|. The disturbed soliton solution 
can be written as 



V; = V;, + = + 0' + = ^n, + 5n e*("^=+'''^) (8) 

where in the linear approximation the real variables (p' and <j)" are related with 
perturbations of the density 5n and the phase 5(j) by the formulae 

, cos(/Ps . „ sinv3„ ^ . 

zon + ^ris smLps ■ dtp, <p — .-on — cos ips ■ dip. (yj 



Substitution of ^ into ([!]) with U(r) ~ and linearizing of the resulting equa- 
tion with respect to (j) give 

i0t + + + - (2|$,p0 + $20*) = 0, (10) 

where asterisk denotes complex conjugation. Let the disturbance be represented 
by a harmonic wave propagating along the soliton: Sn, 64> oc exp[i(p?/ — Q{p)t)] 
OT (j) = [(j)' + icj)") oc exp[i{py — i}{p)t)],(j) = (</>' — #") esc exp[j(py — n{p)t)]. For 
this harmonic form of the disturbance, equation ( 10 ) reduces to 

in^ + iv^c + <(> + h^cc - \p'<t> - (2|<i>.|'0 + $?0*) - 0, (11) 

where C, = x — Vt. Substitution of ([t]) and separation of real and imaginary 
parts yields the eigenvalue problem for i7 = w(p) + iT{p) 



(12) 



ill ii2 w A / r w 
, ^21 i22 y W / V r 

where 

Lii = -VdQ - 2fcT/tanh(fcC), L12 = -5% - 1 + n,, + 2V'^ + 

L21 = i% + 1 - 3n, + 2y2 _ ip2^ L22 = + 2fcytanh(fcC)- ^^'^^ 

It can be solved numerically for any given value of the soliton's velocity V . As 
was found in [7|, the spectrum consist of two parts — unstable = *r(p)) for 
< p < Pc and stable (f2 = ^(p)) for p > p^i where 



Typical plots for r(p) and w(p) are shown in Fig. 1. Non-zero values of F 
for < p < mean instability of dark solitons. The dispersive curve w(p) for 
Pc < p < 00 has an inflection point p = pcr at which 



d uj 
dp'^ 



= (15) 



The found in f5| criterium of transition from absolute to convective instability 
of oblique solitons can be written in the form 
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Figure 1: Growth rate of unstable modes for < p < pc (red) and dispersion law of harmonic 
stable modes for Pc < P (green) as functions of the wave number p of perturbation. Soliton's 
velocity is equal to V = 0.355 and pcr = 1-08. 



where 



Vn 



Mcr COS 6, 



(17) 



is the component of the flow velocity along the soliton, 9 being the angle between 
the direction of the incident flow (x axis) and the oblique soliton, and the 
soliton's velocity V is equal to 



V = Vj_ = Mcr sin t 



(18) 



Hence the dispersion law ui = uj{p, V) depends on Mcr and (16) is actually the 
equation for Mcr as a function of 9. This equation was solved numerically in 
[5] and in the limit of small 9 (deep solitons) the solution tends to the value (|3| 
from below. 



3. Velocity of expansion of unstable disturbances 



Equations ( 15 ) and ( 16 1 can be interpreted in general case, when the soliton's 
velocity V — V± and the flow velocity along it are independent parameters, 
as the statement that the front of the disturbance which breaks the dark soliton 
with formation of vortices propagates into undisturbed region with velocity 
equal to the minimal group velocity 



Vcr = min ( . 



(19) 



If V^i > Vcr, then this wave packet is convected away by the flow from a region 
around any fixed value of the coordinate along the soliton and the instability is 



4 



convective; if V|| < Vcr, then the wave packet evolves into nonUnear stage with 
breaking of the sohton into vortices at any fixed location along the soliton and 
the instability is absolute. 

This formulation can be supported by the following reasoning. At the edges 
of the disturbance its amplitude is small and hence it can be described by the 
linear theory. Therefore we represent evolution of a disturbance as 



(20) 



For large time t the integral can be evaluated by the steepest descent method 
so that the main contribution to it is given by vicinity of the saddle points Ps 
defined as solutions of the equation 



dp 



(21) 



The values of depend ony/t and we suppose that for large enough y equation 
(21) yields two real solutions. Then the integral function includes the fast os- 



cillating factor e'^°^ and the resulting integral corresponds to a usual dispersive 
wave packet slowly decaying with time as However, if with decreasing of 



y the roots of Eq. (21 ) move along the real p-axis and collide with each other at 



some value oiy — yf bifurcating here into two complex roots ps and p* , then the 
saddle points move into complex plane and the integral ( 20 1 includes the factor 
exp[—lm{ps{y))y] exponentially depending on y. Hence, y = y/ corresponds to 
the edge of a large amplitude pulse. At the point of transition of two real roots 



into two complex ones the equation (21 ) has a double root, that is here we have 



d^n 

dp-^ 



0, 



(22) 



which means that the edge of the pulse corresponds to the extremum of the 
group velocity. Just this situation occurs in the case of the dispersion law r2(p) 
of disturbances propagating along a dark soliton, as one can see in Fig. 1, where 
Per corresponds to the minimum ( 19 1 of the group velocity. 



It is instructive to illustrate the above consideration by the example of shal- 
low dark solitons which evolution is described by the Kadomtsev-Petviashvili 
equation [8J and the dispersion law of disturbances was found in explicit form 
by V. E. Zakharov in 0: 



V2 



^{p) = i^V^Jv-Vc, Pc 



V3il-V), 



(23) 



where 1 — V 1. From Eq. (21 ) we find the solution for the saddle point 



2pc + V3^ + ? 



,y 



3^- 2V3p, 



(24) 
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Figure 2: Evolution of a localized disturbance with time. A plane soliton slightly disturbed 
in vicinity of the point y = starts its motion at x = 10 with velocity V = 0.355. The edges 
of a disturbed region expand with velocities equal to Vj = ±1.30. 



which is complex for y^/t'^ < 2pc/VS- Thus, here the disturbance has the 
amphtude proportional to 



0(2/, t) (X exp 



r 

' t 



t2 



Its edge points propagate with velocities 



y± 

t 



±- 



31/4 



±^2{l-V) 



(25) 



(26) 



equal (in absolute value) to the minimum of the group velocity at p = p, 



cr — 3i^c; 





duj 




dp 



(27) 



corresponding to the dispersion law ( 23 1 . 

Thus, we can say that for large enough time the edges of the disturbance of 
a dark soliton propagate into its undisturbed parts with the universal velocity 
Vcr determined by the parameters of the soliton rather than the initial form 
of the disturbance. We have checked this conclusion by numerical solution of 
the GP equation ([T]) with the initial condition in the form of the dark soliton 
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Figure 3: Formation of oblique solitons by the flow of the condensate past an obstacle located 
at the origin of the coordinate system. To avoid generation of the switching wave at t = 0, 
the obstacle potential was turned on gradually during fist 20 units of time. The Mach number 
of the flow is equal to M = 3 and the oblique solitons correspond to the normal component 
of the flow velocity V± = 0.355 equal to velocity of the soliton shown in Fig. 2. The length 
of the oblique solitons grows up with the rate dL/dt = 1 .69 in agreement with the analytical 
prediction dL/dt = 1.71 calculated according to Eq. l|28|. 



Q slightly disturbed at point y = 0. Two stages of evolution are shown in 
Fig. 2. Although the disturbed region evolves fast into the nonlinear stage with 
formation of vortices, its edge points propagate with velocity Vf = 1.30, which 
agrees quite well with the critical velocity Vcr = 1-27 calculated according to 
the linear theory by Eq. ( 19 1 . 



4. Growth rate of the length of obhque solitons 

As was predicted theoretically in [J] for atomic condensate and was confirmed 
experimentally in [5] for polariton condensate, a supersonic flow of condensate 
past an obstacle with size about one healing length generates a pair of oblique 
solitons behind the obstacle provided the incident flow has large enough Mach 
number M > Mcr (see ([s])). If the size of the obstacle is increased, then two 
symmetrically situated fans of dark solitons are generated [ini [H] . After switch- 
ing on the flow velocity (or after the obstacle is put into motion), the length L 
of oblique solitons increases with time linearly and, under supposition that the 
edge of the soliton can be considered as the edge of the disturbance propagating 
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into the undisturbed region, the above theory leads to the conclusion that 

^ = V]i - (28) 

where 

V\\ = M cos e (29) 

is a projection of the flow velocity on the dark soliton direction (hence the 
soliton velocity V is canceled by the normal component V± = M cos 9 of the 
flow velocity and the oblique soliton is a stationary object in the laboratory 



frame of reference). The velocity Vcr is defined by Eq. (19) and represents the 
edge velocity of the disturbance which destroys the soliton into vortices. For 
Vii > Vcr this disturbance is convected away from the obstacle and the oblique 



soliton increases its length L with the constant rate ( 28 1 . 

We have checked this prediction by numerical simulations similar to ones 
performed in [i] but for longer time. Several stages of evolution are shown in 
Fig. 3. As we found, the oblique solitons do not stop their growth even after long 
time of evolution in agrement with the theory [5] and the growth rate agrees very 



well with the analytical prediction (28). In the case of the chosen parameters 



we have V± — 0.355, the component of the flow velocity along the oblique 
soliton is equal to V\\ = Mcos9 = {AP — VlY^'^ = 2.98, and after subtraction 
of the minimal group velocity = 1.27 corresponding to Pcr = 1-08 we get 
dL/dt — 1.71 practically coinciding with the numerical estimate dL/dt = 1.69 
obtained by fitting to the results extracted from Fig. 3. 



5. Conclusion 

In this letter we relate the critical velocity of transition from absolute in- 
stability of dark solitons to their convective instability with the velocity of ex- 
pansion of an unstable wave packet into undisturbed region. A simple estimate 
shows that the expansion velocity is equal to the minimal value of the group ve- 
locity of perturbations propagating along the soliton. This interpretation of the 
absolute/convective instability transition is physically clear and leads to exper- 
imentally verifiable predictions. In particular, the growth rate of the length of 
oblique solitons generated by the flow of BEC past small obstacles can be easily 
calculated. The developed theory can also be applied to description of stability 
of oblique solitons generated by the flow past extended obstacles. For example, 
as was observed in , oblique dispersive shocks generated by the flow of BEC 
past a concave corner start to decay at some moment of time with formation of 
a cloud of vortices. The character of this instability depends on the projection 
of the flow velocity on the inclined side of the corner: if it is greater than the 
minimal group velocity of perturbations, then the instability is absolute and 
the dispersive shock wave will eventually be completely destroyed; otherwise 
the instability is convective and the cloud of vortices is convected away by the 
flow of BEC as long as this cloud does not interact with the boundary of the 
corner and, as a result of this interaction, the convective flow along the corner 
is not destroyed. 
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